Thermalization of nonequilibrium electrons in quantum wires 
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We study the problem of energy relaxation in a one-dimensional electron system. The leading 
thermalization mechanism is due to three-particle collisions. We show that for the case of spinless 
electrons in a single channel quantum wire the corresponding collision integral can be transformed 
into an exactly solvable problem. The latter is known as the Schrodinger equation for a quantum 
particle moving in a Poschl- Teller potential. The spectrum for the resulting eigenvalue problem 
allows for bound-state solutions, which can be identified with the zero modes of the collision integral, 
and a continuum of propagating modes, which are separated by a gap from the bound states. The 
inverse gap gives the time scale at which counterpropagating electrons thermalize. 

PACS numbers: 71.10.Pm, 72.10.-d, 72.15.Lh 
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Introduction. — Relaxation processes and nonequilib- 
rium dynamics in one-dimensional electron systems have 
moved into the focus of recent theoretical [l|-Q and ex- 
perimental research 10l4l3|. One essential characteristic 
of one-dimensional electrons is the absence of relaxation 
in the case of a linear energy dispersion relation and a 
slow relaxation if dispersion is nonlinear, inhibited by 
momentum and energy-conservation laws [4|, 7] ■ Based 
on momentum-resolved tunneling spectroscopy the pe- 
culiarities of relaxation in one-dimesional electron sys - 
tems have been observed in recent experiments 10|, Il2l | . 
Moreover, the violation of the Wiedeman- Franz law ob- 
served [15| at the plateau of the electrical conductance 
in single channel quantum wires can be qualitatively un- 
derstood from the different relaxation processes required 
for equilibration of applied temperature gradients and 
chemical potentials. 

In the presence of many transport channels relaxation 
of electrons at low temperature is primarily provided by 
pair collisions. A solution of the eigenvalue equation of 
the two-particle collision integral in two and three di- 
mensions has already been given four decades ago at the 
early era of Fermi liquid theory [l6| . A remarkable result 
of that theory is that the eigenvalue equation turns out 
to be exactly solvable. In single channel wires, on the 
other hand, conservation laws severely restrict the phase 
space available for scattering processes, and pair colli- 
sions do not provide a relaxation mechanism. If electron 
density is not too low then three-particle collisions [5|, |l7[ 
constitute the leading order relaxation process. Unlike 
the two-particle collision integral in higher dimensional 
systems, a spectral analysis of the corresponding three- 
particle collision integral in single channel wires so far is 
missing. This paper aims to fill this gap. We show that 
for specific momentum configurations of the scattering 
states, relevant for energy relaxation, the collision inte- 
gral of spinless electrons can be diagonalized analytically. 



We find that zero modes of the collision integral are sep- 
arated by a gap from a continuum of decaying modes. 
The gap value provides us the time scale at which coun- 
terpropagating electrons, exposed to a small temperature 
differences, thermalize. 

Formulation of the problem. — Our motivation is to 
identify the time-scale Tth at which thermalization be- 
tween out-of-equilibrium counterpropagating electrons 
occurs in a clean single channel quantum wires. We 
pursue this goal by studying the spectrum of the three- 
particle collision integral under the assumption that left- 
and right-moving electrons inside the wire are initially 
at distinct equilibria, characterized by different temper- 
atures, with AT being the temperature difference. 

Within the Boltzmann kinetic equation approach for 
the electron distribution function, J-' = —I{J-'}, mi- 
croscopic details of relaxation process are stored inside 
the collision integral I{-^}. In the following we spec- 
ify I{J^} for the spinless electrons with quadratic en- 
ergy dispersion, interacting via the Coulomb potential 
V{x) J18i] . In the high density limit interactions are weak, 
e'^/Hvpn <C 1, here vp is the Fermi velocity and k is di- 
electric constant of the host material, and the leading or- 
der relaxation process is due to three-particle collisions: 
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Higher order scattering processes are suppressed by the 
small interaction strength and phase space. In Eq. ([T]) we 
introduced notations, J-i — J-{t, x,pi), and, J-'^ = 1 — J^i. 
The WI2P' = ^\{l'2'3'\VGoV\l23)c\^d{E, - Ef) is 
the rate of three-particle scattering from the incoming 
states Pi,2,3 into the outgoing states Pi',2',3' with en- 
ergies £'i(/) = J2i ^j(j')' respectively. An explicit form 
of W for a generic two-body interaction potential V = 



jLT,k,k,qVgcl^^/^^_^Ck^Ck,, where c^ (cfe) is the elec- 
tron creation (annihilation) operator, has been recently 
derived in Ref. [5|. We here merely mention that Gq 
denotes the free particle Green's function, subscript "c" 
refers to irreducible scattering processes, L is wire length 
and Vq is Fourier transform of interaction potential V{x). 
Because of nonlinearity of the collision integral in 
Eq. Il]) an exact analytical solution of the Boltzmann 
equation is very difficult to find. A simplification is pos- 
sible within the linear response analysis in which distribu- 
tions entering Eq. ([l} can be linearized around an equilib- 
viumstate, Tp = fp+fpf^^p. Here fp= [e(^^-^^)/^+l]-i 
is equilibrium Fermi distribution with ep the Fermi en- 
ergy and temperature T, f^ = 1 — fp, and i/jp in re- 
sponse to externally applied perturbation, in our case 
tpp oc {ep — eF)sgn{p)AT. Restricting to the linear re- 
sponse regime we insert above distribution into the colli- 
sion integral Eq. ([1]) and arrive at the Boltzmann equa- 
tion, ipp-^^ — —C{ip}p-^^ , with the linear collision operator 
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Here the kernel ^{pi};{p.,} 
spectral analysis of the linearized collision operator in 
Eq. ^ under the above formulated assumptions is the 
problem we address in the following. 

Zero modes and symmetries. — Eigenfunctions of the 
collision integral with eigenvalue zero ( "zero modes" ) cor- 
respond to constant in time solutions of the Boltzmann 
equation and are associated with the conserved quanti- 
ties in the system. Indeed, it is readily checked that C 
in Eq. ([2|) is nullified by V'e = £p (energy conservation), 
V'P = V (momentum conservation) , and V'at = const (con- 
servation of total particle number). Since we are inter- 
ested in thermalization we can further restrict to pro- 
cesses in which all of the participating states are close 
to the Fermi points. Then the difference in number of 
left- and right-moving electrons, AiV, is also conserved 
and V'AAf = sign(p) is an additional zero mode. We 
only briefly mention that three-particle collisions chang- 
ing AiV require backscattering and are important for 
the relaxation of differences in the chemical potentials 
of counterpropagating electrons [7[ . 

Defining the Hilbert space of functions endowed with 

scalar product (i/'plV'p) = 2^ I-^ '^P fpfp'^p'^P i* ^^ 
readily checked that C is positive Hermitian, implying a 
spectrum of eigenvalues larger or equal to zero. The zero 
modes form a basis of the four-dimensional subspace of 
conserved quantities. Any tpp that falls off this category 
evolves according to the Boltzmann equation and eventu- 
ally relaxes into one of the zero modes or their linear com- 
bination. In general, the collision operator £ may have 
discrete and continuous parts of the spectrum. However, 
only if the zero-modes are separated by a well-defined 
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FIG. 1: [Color online] Three-particle scattering processes that 
allow for the energy exchange between counterpropagating 
electrons and thus lead to their thermalization. The matrix 
elements (l'2'3'|yG'y|123)c decompose into six contributions, 
A(ll', 22', 33') plus five permutations of primed arguments, 
corresponding to the direct a) and five exchange terms b)-f). 



gap to a smallest nonvanishing eigenvalue the concept of 
a relaxation rate is justified. 

It is helpful to account for an additional symmetry 
of the linearized collision integral. As a consequence of 
invariance of the scattering kernel under reversal of mo- 
menta, A^{pi};{p./} — ^{-Pi};{-Pi/}i ^ commutes with the 
inversion operator H^/jp = 4'-p- Therefore, eigenfunctions 
of C have well-defined parity and the operator itself can 
be decomposed into a direct sum C = £+ © £~ of oper- 
ators C^ acting in the mutually orthogonal subspaces of 
even and odd-parity functions. (Notice that the Hilbert 
space defined above decomposes into a direct sum of even 
and odd-parity functions, mutually orthogonal to each 
other.) As we are interested in the relaxation of an odd- 
parity perturbation, ijj^p = —ipp, our main focus is on 
the spectrum of £^. 

Spectrum, of the linearized collision integral. — We start 
analysis from identifying small parameters in the prob- 
lem. Fermi blocking in combination with conservation 
laws restricts participating scattering states to momen- 
tum strips of order Sp ^ T/vp ^ Pf around the Fermi 
points. Denoting qi — pit — pi the momentum trans- 
fer in a collision we thus have the small parameters 
\(ii/Pj\ ^ T/ep ^ 1, where «,j = 1,2,3. The energy is 
transferred via the three-particle collision in which one of 
the incoming electrons, say with momentum p^, scatters 
off two other counterpropagating electrons (see Fig. [T] for 
the six possible scattering processes where an electron 
with p3 < scatters off two right movers). Using mo- 
mentum conservation 92 + 91 + 93 = we can express 



energy conservation as p2{qi + ga) - (Pi9i + 9? + PsQs + 
9i93 + 93) = which, solved for 53, shows that (73 = 
9i(Pi ^ P2)/2pF + C* ((pi - V2flv\-, ql/py)- It becomes 
clear now that energy transfer between the counterprop- 
agating electrons occurs via small portions of momen- 
tum exchange VF\q3\ ^ T'^/ep, and {\q3,/qi\, \q'i/q2\} < 1 
present additional small parameters, while qi ~ —52 up 
to corrections of O (92(^1 — P2)/pf))- 

A further important ingredient for our calculation is 
the three-particle scattering rate for the Coulomb poten- 
tial. Expanding Wl^^^' = K^^^^' 8{E, - Ef)6g,+g^+g^=o 
to the leading order in the above small parameters, we 
find that the scattering rate depends only logarithmi- 
cally on the momenta transfer qi. Indeed, Aj23"^ = 

x(¥)^[t^]'r2('zi-'?3),with 



extend solutions to negative values p < by taking even 
and odd parity combinations V'p = 0{p)'4'p ± 0{—p)'ip-p- 
We next show that Cq can be transformed into a lin- 
ear second-order differential equation. As a first step, 
we adopt logarithmic accuracy approximation to substi- 
tute the argument of the logarithm in the momentum 
dependent scattering rate Eq. ^ by its typical value, 
q\/{qsPF) — > 1, as dictated by the conservation laws. The 
linear collision operator with constant scattering rate can 
be reduced then to 
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r(gi,g3) = [l + i{lE + HkFw])]\n[2pF\q:i\/ql]. (3) 

where 7^ is the Euler constant. This result [Eq. ([3])] ap- 
plies at not too low temperatures {Td/hvp, kpd} ^ 1 
at which screening of the nearby gate is not impor- 
tant. Technically speaking, the typical Fourier compo- 
nent Vfc of the interaction potential then has wave num- 
ber in the range d~^ ^ |fc| ^ w~^, so that Vk — 
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-Jk\^ ^ -^r- -^ -^ifr\ M- The bgarithmic q,- 
dependence of the scattering rate is a result of subtle 
cancellations between direct and various exchange terms 
contributing to the scattering amplitude, see also Fig. [T] 
It is owing to this property of W that makes an analytical 
diagonalization of C possible. 

Our strategy is now to split the linearized collision op- 
erator into two contributions 



cm = Com + scm , 



(4) 



where Co allows for an analytical diagonalization, while 
corrections to the spectrum from SC turn out to be 
irrelevant for our problem (and may, in principle, be 
calculated from perturbation theory). More precisely, 
we choose £0 as in Eq. 1^ with kernel /Cq resulting from 
the original /C upon linearization of the quadratic energy 
dispersion, i.e. upon substituting fp±pj, and S{Ei — Ef), 
respectively, by.g±p = (e±^^P/^+l)-i, emdSoiE.-Ef) = 
^itr; ['5gi-l-g2^g3 0(Pi>P2, -Pz) + <5g2 +93^916(^1. "^2, ~Pz% 

e(pi,P2,-p3) = \Q{PX)B{P2)Q{-PZ)+Q{-P1)K-P2)Q{PZ)\- 

That is the operator £0 corresponds to a problem with 
linear dispersion relation. The expectation that correc- 
tions to the spectrum from bC are small is, of course, 
a consequence of the fact that three-particle scattering 
processes which provide energy exchange involve all col- 
liding electrons near the Fermi points, while states at the 
band bottom are not crucial for thermalization. Since 
the above approximation preserves inversion symmetry, 
[^OiII] = 0, eigenfunctions of Co have well-defined par- 
ity. Given this property we may restrict the eigenvalue 
equation i3o{V'"}p = '^nV'p to momenta p > and then 



where k = i^T jvp is characteristic momentum due to 
thermal smearing of Fermi functions, hp = (e'"^^/"^ — 1)"^ 
is the bosonic distribution function, g^ = 1 — 5^, and 70 — 

iJ^TL^^2^Y\^j^nAkF^\- The sequence of 
exact transformations leading to Eq. ([5]) can be sum- 
marized as follows [18|. (i) It is convenient to orga- 
nize Co into six contributions Co = J2s=± S ;= 






where l"- iih} = ^-7; 
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/Cq [4ij — V'j'] and 

1Cq~ {1Cq~) describes processes in which the right- mover 
Pi scatters off one right- and one left-mover (two left- 
movers), (ii) In the individual contributions H" we re- 
move two out of five momentum-integrations employ- 
ing conservation laws and complete two further integra- 
tions with help of the identities: X]i)5p.9io-g — hl^i^ 

yq9p+q9p-q = -j^Pb-p. and yqqgp+qhq = -^(^^ + 

p^)gp- (iii) We find that all terms 17 ~ {i^} and l^ {-0} 
are identical zero, while I'l {0}+^^ {V'} can be summed 
to give Eq. (O . 

We then introduce odd and even momentum combina- 
tions with respect to the Fermi point 



^P = \l9p9p [V'pF+p ± i^PB 



(6) 



where the normalization factor is chosen for convenience, 
and recast the eigenvalue problem for odd and even com- 
binations, s = ±, as follows 
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with 6± — 0, 1, respectively. At this stage we intro- 
duce dimensionless momentum k = ^^^, energy r2„ — 

-2^ (t^) ^n, and notice that the kernel in Eq. ^ de- 
pends on the difference of its arguments, which makes it 
convenient to perform a Fourier transformation ip' 
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cosh'^ X ' 



As a result, the 
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eigenvalue equations for the Fourier images "0^ * reduce to 
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Schrodinger equations of a particle moving in a Poschl- 
Teller potential 
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which can be solved with help of operator-algebra tech- 
nique known from the harmonic oscillator problem [19] . 

The eigenvalue problem in Eq. (jS]) allows for one even- 
and one odd-parity bound-state Qn = of the form 
■!/)°+ = 1/ cosh a: and -0°" = — 3sinha;/cosh^ x, respec- 
tively |2C| . Upon inverse Fourier transformation and ex- 
tension to negative momenta in the above prescribed 
manner this gives the four zero-modes ipj^/ = const, 
^AN = sgn(p), ipp^ = p, and ip°j^ = \p\. As already 
discussed, the first three functions are consequences of 
conservation of total particle number, momentum and 
the difference in number of left- and right-moving elec- 
trons. The fourth zero mode expresses conservation of 
energy for the linearized spectrum. Of course, these four 
zero modes could have been directly inferred from Cg. 

More relevant for our problem is the fact that these 
bound-state solutions are separated by a gap Sfl — 1 
from a continuum of propagating modes and Co, there- 
fore, possess a well-defined smallest nonvanishing eigen- 
value. In order to associate this latter with the ther- 
malization rate, we have to make sure that this gap is 
also present in the original collision operator C~ . Em- 
ploying that eigenfunctions of Co form a complete set we 
may express C~ in terms of its odd-parity subset. Fi- 
nally, reminding that tpAN and ipp are also nullified by 
C~ , it is evident that the smallest nonvanishing eigen- 
value of £^ is of order Slo = min ,- {(-0^1/3 iV'jT/)}, with 
{tpjl } the set of eigenstates corresponding to odd-parity 
propagating solutions of Eq. (|8]) . Since matrix elements 
\{i:^\SC\^^,)\ < \{^k\^o\^k')\ it readily follows that C" 
and Cq share a gap of same order. Employing then 
Sfl — 1, restoring original units, and inserting the ex- 
plicit form of 7o we arrive at the thermalization rate 



1/rth = c{eF/K){e^/hvFnf\^{kFw){T/eFf 



(9) 



where coefficient c = In 2/27r^ and function \{x) = 
x"^ \a.x. Equation @ presents the main result of this pa- 
per. We note here that the temperature dependence of 
T^ can be understood from a simple phase-space argu- 
ment like in Fermi liquid theory. To first approximation 
the thermalization rate follows from the out-scattering 
part of the collision integral, l/rout oc X]|n > ^{pi}' where 
two out of the five momenta {pi\ are fixed by the con- 
servation laws, while the remaining three extend over 
momentum range set by the temperature broadening 
5p ~ T/vf of Fermi distributions. Since the scatter- 
ing amplitude is only logarithmically dependent on mo- 
mentum transfer then scattering rate W ~ const and, 
therefore, l/rout oc {T/sFf ^■ 



Conclusions. — We have analyzed the spectrum of the 
three-particle collision integral in a one-dimensional elec- 
tron system. We found that zero modes, associated with 
the conservation laws, are separated by an energy gap 
from a continuum of propagating modes, and identified 
the gap with the relaxation rate rth relevant for ther- 
malization of counterpropagating electrons. Our analysis 
applies to clean single channel quantum wires of spinless 
electrons at not too low densities. It is highly desirable, 
yet very challenging, to extend the analysis to lower den- 
sities where interactions become strong IJ]. Also, inclu- 
sion of spin [j| presents an interesting problem, since at 
very low temperatures the effect of spin-charge coupling 
becomes relevant [23]. 
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